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Stratified disks are locally stable 
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Abstract. Notwithstanding recent claims by Richard et al., there is no linear hydrodynamic instability of axisym- 
metrically stable disks in the local limit. We prove this by means of an exact stability analysis of an unbounded 
incompressible flow having constant stratification and constant shear. 
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1. Introduction 



In a very recent preprint, Richard et al. (2001 ) claim to 
demonstrate nonaxisymmetric local incompressible linear 
instability of shearing disks with vertical stratification. 
(They also discuss more general flows that are strati- 
fied in both radius and height). The growth rate is ap- 
proximately the geometric mean of the shear rate (S) 
and Brunt- Vaisala, frequency (N), hence quite rapid since 
| S/Q\ = 3/2 in keplerian disks, and typically one expects 
N ~ fl away from the disk midplane. If true, this would 
be very important, as the instability might support tur- 
bulence and angular-momentum transport in disks too 
weakly ionized for magnetorotational instabilities ( Balbus 
|fc Hawley, lggf ). 

In we derive the exact equation ( |l4| ) for time- 
evolution of spatial Fourier modes in a Boussinesq flow 
with a linear shear profile. In we use this equation to 
show that all local disturbances are bounded at late times, 
and hence there is no growing local mode. We also briefly 
discuss how global instabilities might arise in the presence 
of reflecting boundaries. 



2. Local Boussinesq model 



Using a local WKB analysis, Richard et al. (2001) assert 
that local Fourier modes with wavenumbcr k = [k r , kg,k z ) 
grow oc exp(rt) [j] 



r z = _ 



2klN 2 ns 



k 2 {N 2 + 2QS) + k 2 N 2 + k 2 



(1) 



1 The term k^N 2 in the denominato r is missing in the ex- 
pression quoted in Richard et al. (2001 ). This appears to be a 
typographical error. Whether the term is present or not does 
not affect our critique. 



As usual, ril(r)ee is the background orbital velocity, S 
dfl/dlnr < is the radial shear, and 



,d{r 4 n 2 ) 
dr 



20(2f2 + S), 



J_<Wd_ /P_ 

■fp dz dz \p~< 



are the squares of the epicyclic and Brunt- 
Vaisala frequencies. Since it is assumed that N 2 > 
and k 2 > 0, the flow would be stable against both 
convection and centrifugal instability if there were no 
shear. Yet according to eq. ([!]), there is a constant growth 
rate to ~ y— 2Q.S in the limit of strong stratification 
A^ 2 — * oo provided only that d(fl 2 )/dr < (as is normally 
the case in astrophysical disks). 

In the extreme local limit, one considers disturbances 
on scales very small compared to the distance to the near- 
est boundary of the flow. Also, provided that f2 and A^ 2 
are smooth functions, they vary little on the small length- 
scale. It is therefore reasonable to consider an unbounded 
flow in which all of the frequencies £1,N,S,k are constant 
in space and time. Variations in density can be ignored 
except in the buoyancy term (Boussinesq approximation) . 
The governing equations for this flow in a corotating frame 
are 





L 




V-v. (2) 



The quantity $ represents the fractional density difference 
between a displaced fluid element and the background 
(which the Boussinesq approximation considers to be in- 
finitesimal) . L is a constant scale height for the entropy, so 
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that N 2 = g/L if g is the vertical gravitational accelera- 
tion, if) is a scalar potential for accelerations due to gravity 
and pressure. The background state is 

v Q = Sxe y , i? = 0, = QSx 2 . (3) 

The coordinates r, 6*, z have been replaced by local 
Cartesians x,y,z. Were Richard et al. (2001)'s analysis 
correct, their growth rate (|l|) would apply to this model. 

At any given time (t) , linearized disturbances of this 
flow can be be decomposed into spatial Fourier compo- 
nents, 



vi = eVe 1 



0i = e0e* 



0i = e^e 1 



(4) 



where e < 1 is a formal small parameter. Following a 
method attributed to Kelvin (1887), (^) is a solution to 
the linearized equations of motion for all time if V, 0, \P 
emd fc are appropriate functions of t [but independent of 
position r = {x, y, z)]. Note that this means that the time 
dependence of the flow attributes will not, in general, be a 



simple exponential function, in contrast to the Richard et 
al. (20([l|) assumption. Substituting into the first of eqs. (|2|) 
and collecting terms of first order in e, 







dk 
~dt ' 
dV 
~dt 



Sxk v iV 



SV x e v + 2Vle z x V + ik^> + iV L6e. 



The coefficients of each coordinate xyz must vanish sepa- 
rately, whence 



S 



dt 



SV x e y + 2tte z x V + ik^ + N 2 L&e z 



~dT 

and 
dV 
~dt 

By similar steps, 
d& _ V, 
~dt ~ ~L 
and 

k V = ^ V z = - 



dk z 
~~dt 



0. 



0. 



(5) 



(6) 



(7) 



k x V x 



kyVy 



Combining the time derivative of 
* = A [2(5 + n)k y V x - 2VLk x V y + N 2 Lk z Q] 



(8) 

) with eq. (||) yields 



Note the abbreviations 



(9) 



(10) 



Since V z can be eliminated via eq. (jg), the system (g)- 
(0) is third order in time.^J However, it has a first integral, 



kyU x 



(20 -I- S)k z LQ = constant 



(11) 

Yet Richard et al. (2001 )'s dispersion relation is fourth- 



order in frequency. At k y — 0, two roots vanish, and the claim 
is made that these develop into ±r as given by ([!]) when k y ^ 0. 
Non-vortical perturbations, however, are restricted to obey a 
second order system. Moreover, if time dependences were in 
fact exponential, then a nonzero value of the vorticity constant 
in ([ll]) would correspond to a zero eigenfrequency at any k y . 



representing conservation of potential vorticity. For grow- 
ing modes, the constant would be asymptotically negligi- 
ble. So we set it to zero. Eliminating V z between eqs. ^ 
and (§): 



k x U x + kyUy — k z L——. 



dt 



We use eqs. ([U| & (|12j 
eq. (||) in favor of 9 & 

d 2 B 
dt 2 



= 



2 S kg>ky kg 



k 2 

T 2 _L 



k 2 k\ 



(12) 

to eliminate U x and U y from 

de 

~dt 



N 2 ^ + n 2 -^ + 2S(2Q- 
k l k 



S) 



k 2 k 2 
k 2 k\ 



0. 



(13) 



(Note that eq. ( |l3| ) also describes the evolution of vertical 
Lagrangian displacements, a consequence of eq. (0).) The 
variable 

$ = A e 

obeys a slightly tidier equation: 
d 2 $ 







d< 2 



k 2 k 2 
k 2 k 2 ( h 2 




(14) 



Together with the solution of eq. (^), 

k x (t) = k x (Q) — 2Stk y , k y ,k z — constants, (15) 

eq. ( |l4| ) governs the evolution of a Fourier mode for the 
density perturbation. The equation is exact within the 
context of the model defined by eqs. (|^)-(^). In fact, it 
is not even restricted to small amplitudes, because the 
terms quadratic in e vanish when eq. (^) is substituted 
into cq. (||) as a consequence of eq. (||) . 

3. Discussion 

Clearly, eq. (|l|) is incompatible with eq. (fl4|). When k z 
0. for example, eq. (Q) predicts a finite growth rate. The 
exact solution of eq. (|lj) corresponds to vertical oscilla- 
tions at the Brunt- Vaisala frequency, an elementary result 
(consider purely vertical displacements) that may be ob- 
tained directly from the z component of equation (JsJ) - 

More generally, as St — > ±oo, fc 2 « k 2 w k 2 > k 2 ,k 2 , 
so that 0J 2 (t) — > N 2 and $ cx exp(±zAt). Hence , contrary 
to the result claimed by Richard et al. (2001 ) [eq. (0) 
there is no exponentially growing linear modei 



3 This assumes that the disk is stably stratified. On the 
contrary, when N < 0, an exponentially growing convective 
disturbance persists as k x —* oo if, as here, viscosity and t her- 
mal conduction are ignored {e.g., Ryu & Goodman, 1992). 
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An approximate dispersion relation can be read off 
from eq. (14) in the limit k x 3> k y ,k z (i.e. at late times). 



The corresponding radial group velocity is 

fa ^ S(N 2 ~K 2 ) kjky 

dh x ~ + N kl [ 1 

This tends to zero as t~ 3 . It can be shown that a radi- 
ally localized wavepacket, obtained by superposing Fourier 
modes with a range of k x (0), asymptotically approaches 
resonant radii where the frequency of the disturbance 
= ±N in a frame corotating with the local background. If 
the disturbance is excited with azimuthal wavenumber k y 
by a force corotating steadily at radius tq , the resonances 
lie approximately at ro ± N/Sk y . In order that our local 
analysis apply to radially bounded flows, it is necessary 
that these resonant radii lie inside the fluid. Otherwise 
the wavepacket may reflect from the boundary, suffering 
a change in the sign of k x , and offering the possibility of 
repeated passages through the "swing amplifier" at small 
k x /k y where to 2 can be briefly negative.^ 



Th e instability discuss ed by by Molemaker 
McWil[ iams, fe Yavneh (200l| ) for s tratified, centrifugally 



stable Couette flow, and cited by Richard et al. (2001 ) 
as support for their claims of instability, is an intrinsi- 
cally global phenomenon, and has no local counterpart. 



Molemaker, McWilliams, fc Yavneh (2001 ) analyzed their 



instability in terms of trapped Kelvin wave edge modes; 
it is also possible that the resonant interaction described 
above may be involved in destabilizing some of the modes. 
In either case, the mechanism is certainly not reducible 
to exponentially-growing local Brunt- Vaisala, oscillations. 
A full analysis of the global problem is beyond the scope 
of this note. 
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See [Toomre ( 198l] ) and |Goldreich fc Tremaine (1978|) for 



particularly clear discussions of the swing-amplifier in the con- 
text of two-dimensional compressible disks. 



